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COMMENTS  ON  CERTAIN  BOUNDARY  LAYER  PARAMETERIZATION 
SCHEMES  USED  IN  ATMOSPHERIC  CIRCULATION  MODELS 

1.  INTRODUCTION 

The  structure  of  the  atmospheric  boundary  layer  is  characterized  by 
two  distinct  sublayers.  Immediately  above  the  surface  up  to  approximately 
30  m height  the  vertical  variation  of  eddy  stress  and  other  fluxes  can 
be  neglected,  and  the  wind  direction  is  essentially  constant.  This 
layer  is  called  the  surface  layer,  or  alternatively  the  constant  flux 
layer.  The  height  of  the  constant  flux  layer  is  quite  variable.  At 
night,  when  fluxes  are  weak  and  the  stratification  strong,  the  stress 
may  be  constant  for  only  a few  meters  near  the  surface.  At  the  other 
extreme,  during  the  daytime  convective  periods,  there  is  strong  tur- 
bulent mixing.  Thus  the  surface  layer  may  extend  to  a considerable 
height  above  the  surface.  Above  the  constant  flux  layer  is  the  Ekman 
layer,  or  the  outer  layer.  In  this  layer,  the  earth's  rotation  and 
the  large  scale  pressure  gradients  in  addition  to  the  eddy  stresses 
are  important  in  controlling  the  wind  and  temperature  profiles.  The 
turbulent  fluxes  typically  decrease  upward  within  the  outer  layer. 

The  interaction  between  the  boundary  layer  and  the  free  atmosphere 
occurs  through  the  vertical  fluxes  of  momentum,  mass,  heat  and  moisture, 
and  the  frictionai  dissipation  within  the  boundary  layer.  The  essential 
purpose  of  boundary  layer  parameterization  is  to  generate  simple 
models  for  the  computation  of  the  boundary  layer  fluxes.  The  first 
part  of  this  paper  discusses  the  parameter i zat ions  based  on  the  surface 
layer  similarity  theory.  This  theory  is  applicable  whenever  the  first 

Note:  Manuscript  submitted  June  23,  1977. 
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grid  level  of  a general  circulation  model  can  be  placed  within  the 
constant  flux  (or  surface)  layer.  In  particular,  we  shall  examine  the 
effect  of  varying  the  constant  flux  layer  height  on  the  computation  of 
the  boundary  layer  fluxes.  The  second  part  of  this  paper  concerns 
parameterization  of  the  whole  planetary  boundary  layer.  This  is 
neccess i tated  by  the  fact  that  many  large  scale  models  have  such  poor 
vertical  resolution  that  the  first  interior  grid  point  level  is  well 
above  the  top  of  the  surface  layer.  The  boundary  layer  fluxes  in  this 
case  are  calculated  by  the  so-called  Rossby  number  similarity  theory  or 
generalized  similarity  theory  with  empirically  determined  similarity 
functions.  According  to  the  similarity  theories,  these  empirical 
functions  depend  on  the  surface  Mon  in-Obukhov  length  which  is  impractical 
to  determine.  However,  a bulk  Richardson  number  which  is  closely  re- 
lated to  the  Mon i n-Obukhov  length  can  be  explicitly  determined.  Thus, 
we  shall  specifically  discuss  how  the  similarity  functions  can  be 
determined  by  the  bulk  Richardson  number  approach.  Moreover,  we  shall 
compare  the  computed  boundary  layer  fluxes  of  the  generalized  similarity 
theory  and  the  Rossby  similarity  theory. 

2.  PARAMETERIZATION  BASED  ON  THE  SURFACE  LAYER  SIMILARITY  THEORY 


A.  Constant  F I ux  Layer  Theory 

In  this  section,  we  shall  discuss  computation  of  surface  fluxes 
of  heat  and  momentum  based  on  flux-profile  relationships  as  established 
through  the  Mon i n-Obukhov  similarity  theory.  The  most  complete  char- 
acterization of  fluxes  and  profiles  comes  from  the  Kansas  experiments, 
as  published  by  Businger  et  al.  (1971)- 

Following  Businger  et  al . (1971),  we  define  the  following  dimension- 


less  quantities: 
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^ q80/3Z 

Richardson  number,  a 

' 0(3u/9Z)^ 

stability  parameter 

r 

* kZ  3u 

"^m  “ u,,  3Z 

A dimensionless  wind  shear 

k 

A dimensionless  temperature 

^0  0^  3Z 

grad  lent 

kqw‘e'2 


A dimensionless  height 


where 

g = 

k = 

U,  V,  w = 
u = 

Z = 

0 » 

U...  = 

L = 

P ” 

"T  S 

‘o 

c « 

p 


acceleration  due  to  gravity 
von  Karmans  constant  = (0.35) 

longitudinal,  lateral  and  vertical  components  of  the  wind 

magnitude  of  the  mean  horizontal  wind  vector 

vertical  coordinate 

potential  temperature 

friction  velocity  {(x^/p)^} 

scaling  temperature  {-w' 6 ' / (ku^) } 

Obukhov  length  {-9u.?,/ (kgw' 6 ' ) } 
air  density 

surface  shearing  stress 

specific  heat  of  air  at  constant  pressure 


Based  on  the  experimental  results  of  Businger  et  al.  (1971),  the 
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universal  functions  d)  and  (j)  can  be  written  as; 

m o 


(1  - 150'^  . 

' < 0 

(unstable) 

I + 

C > 0 

(stab  1 e) 

0.7^(1  - 90'^, 

A 

0 

(unstable) 

0.74  + 4.75  , 

? > 0 

(stable) 

The  procedure  involves  the  computation  of  the  surface  shearing  stress 

T and  heat  flux,  pC  Uj.9 , , for  a given  set  of  wind  and  temperature  pro- 
o p -- 

files  between  Z and  surface  roughness  height  Z^.  The  height  of  the 
surface  or  constant  flux  layer  Z is  quite  variable  as  it  was  pointed 
out  previously.  As  a result,  an  assumed  height  of  the  surface  layer 
can  have  great  effects  on  the  computation  of  the  surface  shearing  stress 
and  heat  fluxes.  Here  we  attempt  to  assess  the  effects  on  the  surface 
fluxes  of  momentum  and  heat  based  on  experimental  data  when  different 
constant  flux  layer  heights  are  assumed. 

B.  Analysi s Procedure 

Data  used  for  the  analyses  are  based  on  the  Wangara  (Clarke 
et  al.,  1971),  and  the  Kansas  (Izumi,  1968)  field  experiments.  From 
the  Wangara  data,  AO  days  of  wind  and  temperature  measurements  at  four 
levels,  i.e.  50  m,  4 m,  2 m,  and  1.2  m were  chosen.  From  the  Kansas 
data,  hourly  values  of  wind  and  temperature  at  five  levels,  i.e.  32  m, 

16  m,  8 m,  4 m,  and  I m were  selected. 

Only  a brief  account  of  the  computation  procedure  of  the  surface 
fluxes  will  be  given  here.  For  details  the  reader  is  referred  to 
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Barker  and  Baxter  (1975).  First,  the  bulk  Richardson  number  was  cal- 
culated from  the  wind  and  temperature  profile  data.  Next,  the  surface 
Obukhov  length  was  deduced  on  the  assumption  that  an  individual  ver- 
tical flux  and  its  associated  vertical  gradient  strictly  obey  the 
empirical  relationships  of  Businger  et  al . (1971).  Then,  the  surface 
friction  velocity,  and  thus  surface  shearing  stress,  were  calculated 
from  the  integrated  form  of  the  empirical  diabatic  wind  shear  relation- 
ship, using  observed  wind  speeds  at  selected  heights.  It  should  be 
noted  that  no  assumption  was  made  on  the  value  of  the  surface  roughness 
length.  La  the  heat  flux,  and  thus  the  kinematic  heat  flux,  are 

obtained  lux  gradient  relationship  of  Businger  et  al.  (1971) 

using  th  calculated  values  of  friction  velocity. 

C.  Analysis  Resul ts 

J 

Fig.  1 shows  the  results  of  the  analyses  based  on  the  Wangara  I 

data.  Eight  day's  results  are  presented.  These  are  typical  of  the  40 
day's  results  of  the  Wangara  data.  The  most  striking  features  occur 
during  the  convective  hours,  i.e.,  1200  and  1500.  The  shearing  stress 
(top)  and  heat  fluxes  (bottom)  calculated  using  wind  and  temperature 
profiles  at  Zj  = 1.2  m and  = 50  m show  great  discrepancies  when 
compared  to  those  calculated  with  profiles  taken  from  Zj  * 1.2  m and 
Z2  “ 2 m.  We  speculate  that  if  wind  and  temperature  profiles  at  a 
height  of  50  m (or  greater)  were  used,  the  very  super-ad iabat ic  layer 
existing  near  the  ground  during  the  convective  hours  may  not  be  re- 
flected in  the  calculation  of  a bulk  Richardson  number.  As  a consequence, 
the  Obukhov  lengths  are  calculated  to  be  of  nearly  neutral  values,  and 
thus  the  fluxes  calculated  are  much  too  small  in  magnitude.  However, 
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if  wind  and  temperature  at  a height  lower  than  2 m are  used,  then 
measurement  errors  in  the  data  will  probably  lead  to  inaccurate  flux 
calculations.  Note  that  both  shear  stress  and  heat  fluxes  calculated 
with  Zj  = 1.2  m and  ^ m compare  well  with  those  of  Melgarejo  and 

Deardorff  (1975). 

The  1968  Kansas  surface  layer  data  include  independent  measurements 
of  both  surface  friction  velocity  and  heat  fluxes.  Thus,  the  calculated 
results  can  be  easily  compared  with  the  measured  values.  Four  cases  are 
considered.  Case  A:  the  height  of  the  constant  flux  layer  (CFL)  is  set 
to  be  A m with  wind  and  temperature  profiles  taken  at  two  levels  i.e., 

Z,  = 1 m,  and  Z_  = ^ m.  Case  B:  CFL  = 3 m,  Z,  = 1 and  Z.  = 8 m. 

I Z I m Z 

Case  C:  CFL  = 16  m,  Z,  = 1 , Z.  = 16  m.  Case  D:  CFL  = 32  m,  Z,  = 1 m 

I m Z I 

and  Z^  = 32  m.  The  calculation  of  surface  friction  velocity  and  kinematic 
heat  flux  are  done  in  the  same  way  as  for  the  Wangara  data  described 
previously. 

The  calculated  values  for  these  four  cases  for  the  Kansas  experiment 
are  tabulated  in  Appendix  I.  The  RMSE  values  (Root  Mean  Squares  Errors) 
of  the  calculated  surface  friction  velocity  and  kinematic  heat  flux  with 
respect  to  the  measured  values  are  given  in  Table  1.  It  can  be  seen 
from  a comparison  of  the  RMSE  values  given  in  Table  1 that  the  poorest 
results  were  obtained  when  wind  and  temperature  profiles  were  taken  with 
Zj  “ 1 m and  Z2  = ^ m (Case  A).  On  the  other  hand,  the  best  results  in 
the  calculation  of  surface  friction  velocity  and  kinematic  heat  flux 
come  from  Case  D where  wind  and  temperature  data  at  Z^  = 1 m and  = 

32  m were  taken.  These  results  again  suggest  that  if  the  height  of  a 
constant  flux  layer  is  assumed  to  be  very  shallow  such  as  less  than  k m 


i 

I 


6 


depth,  then  any  small  measurement  errors  in  wind  and  temperature  profiles 
will  lead  to  great  inaccuracies  in  the  calculation  of  surface  shearing 
stress  and  heat  fluxes. 

The  results  of  the  above  analyses  are  oaseJ  on  an  assumption  that 
the  data  taken  from  the  Wangara  and  Kansas  field  experiments  strictly 
followed  the  flux-profile  relationships  of  Businger  et  al . (1971).  A 
more  stringent  assumption  was  that  all  the  data  taken  at  any  two  levels 
would  behave  according  to  the  surface  similarity  theory.  These  two 
assumptions  might  be  responsible  for  the  fact  that  great  variability 
was  exhibited  in  the  surface  flux  calculations  with  different  constant 
flux  layer  heights  assumed  in  the  analyses.  For  the  purpose  of  para- 
meterizing the  atmospheric  boundary  layer  for  use  in  general  circulation 
models,  a fixed  height  of  10-30  m for  the  constant  flux  layer  may  be 
practical  and  reliable.  However,  for  modeling  detailed  structures  of 
the  atmospheric  boundary  layer,  a sensitivity  analysis  of  the  effect 
of  different  CFL  heights  is  a necessity. 

3.  PARAMETERIZATION  BASED  ON  THE  MATCHING  OF  SURFACE  AND  OUTER 

LAYER  SIMILARITY  THEORIES  ! 

As  pointed  out  previously,  when  a GCM  has  such  a poor  vertical 
resolution  that  the  first  interior  grid  level  is  well  above  the  top 
of  the  surface  layer,  one  must  parameterize  the  whole  boundary  layer. 

The  best  forms  of  parametric  relationship  for  this  have  been  derived 
by  the  matching  of  mean  profiles  predicted  by  surface  and  outer  layer 
similarity  theories  (Blackadar  and  Tennekes,  1968;  Z i 1 i t i nkev i ch , 1975). 

The  general  form  of  the  drag  (or  surface  shear  stress)  and  heat  transfer 
(or  heat  fluxes)  relations  obtained  from  similarity  matching  arguments 

ui 
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is, 

ku./u,  = - (2,nZ  + A) 

n " o 

kv^/u.,.  = - B sign«f 

k(e,  - 9„)/6,  = - (£nZ  + C) 
n Q " o 

k(q|^  - D) 

where  u and  v are  the  horizontal  components  (in  the  direction  of  surf-  e 
shear  and  perpendicular  to  it,  respectively)  of  mean  velocity  vector 
(the  subscript  h refers  to  the  variables  at  the  top  of  the  boundary 
layer),  Z^  is  the  roughness  parameter  normalized  by  the  scale  height 

(h  or  u.,yf)  of  the  boundary  layer,  and  A,  B,  and  C and  D are  some 

similarity  functions. 

A.  Determinat ion  of  the  Similarity  Funct i ons 

The  purpose  of  this  section  is  to  describe  how  the  similarity 
functions  A,  B and  C can  be  determined  by  the  bulk  Richardson  number. 

In  theory,  the  similarity  functions  A,  B and  C are  functions  of  at- 

mospheric stability  u,  where  h is  height  of  the  planetary  boundary  • 

2 

layer  and  L is  the  surface  Mon  in-Obukhov  length  defined  as,  L = 6u.,.  /kge 
It  is  apparent  that  L is  related  to  the  internal  parameters  u.,.  and  0_,. 
which  are  to  be  determined.  Thus,  it  is  not  practical  to  determine 
similarity  functions  A,  B and  C in  terms  of  L.  Deardorff  (1972)  dis- 
cussed the  relationship  between  the  geostrophic  drag  or  heat  transfer 
coefficients  and  the  bulk  stability  (bulk  Richardson  number).  Recently, 
Yamada  (1976)  reexamined  the  similarity  functions  based  on  the  Wangara 
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data.  Both  in  the  stable  and  unstable  conditions,  his  proposed 
universal  functions  A(u) , B(u)  and  C(u)  show  substantial  improvement 
when  compared  with  the  observed  data  over  the  previous  work  by 
Deardorff  and  Melgarejo  (1975)  and  Arya  (197^).  We  shall  thus  re- 
formulate the  similarity  functions  in  terms  of  bulk  Richardson  number 
based  on  Yamada's  (1976)  work. 

The  bulk  Richardson  number  R.„  is  related  to  similarity  functions 

I D 

by 


^ Hn(h/Z^)  - C 

*^iB  " r {in(h/Z  ) - A^}  + 

o 


where  functions  A,  B and  C are  defined  as  in  Yamada  (1976).  It  follows 
that  for  a given  set  of  values  of  ? = h/L  and  roughness  parameter  Z^, 
the  corresponding  values  of  A,  B,  and  C,  and  R.  are  readily  determined. 

I D 

Now,  for  a given  values  of  R.„  and  Z , we  like  to  determine  h/L  and 
’ I B o 

hence  A,  B,  and  C.  For  this  purpose,  the  Newton-Ralphson  method  is 
usefu 1 , 


where 


and 


f(^)  = 0.74r,(a  - c)/{(a  - A)^  + B^} 


f(<;)  = f(;.){i/;  - c'/(a  - c)  - z(aa'  - «a'  + bb')/ 


(a  - A)^  + B^  } 
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where  <;  = h/L  and  £n(h/Z^)  = a.  The  subscript  i denotes  the  interation 
index,  and  { )'  the  derivates  with  respect  to 

The  Newton-Ra 1 phson  method  was  tested  both  for  stable  and  unstable 
cases  of  similarity  functions  A,  B,  and  C based  on  Yamada  (1976).  We 
found  for  stable  cases,  i.e.  h/L  > o,  the  Newton-Ra 1 phson  method 
converges  within  five  iterations.  Following  Yamada  (1976),  the  sim- 
ilarity functions  for  stable  cases  are: 

(i)  18  > ? > 0 


A = 1.855  - .380^,  A'  = - .380 

B = 3.020  + .300c,  B'  = .300 

C = 3.665  - .819c,  C = - .819 


(ii)  35  > C > 18 


A = 1.855  - .380c,  A'  = - .380 
B = 3.020  + .300c,  B'  = .300 
C = - i*.32(c  - 11.21)^,  C = - 4.32(c  - 11.21)'^ 

(iii)  ; > 35 

A = - 2.9^(C  - 19.9^)^,  A'  = 1.47(c  - 19.9^)'^ 

B = 2.85(C  - 12.47)^,  B'  = 1.425(C  - 12.47)'^ 

C = - 4.32(c  - 11.21)^,  C = 2.16(C  - 11. 21)"^ 

For  unstable  cases,  a simple  linear  relationship  occurs  between 
the  bulk  Richardson  number  and  C.  i.e.. 


5-  f(C,) 


where  C„  = t £n(h/Z  ),  where  k is  von  Karman  constant,  and  f(C,,)  = .75 
N k o N 


C, , - 1.  Table  2 shows  the  values  of  A,  B,  and  C based  on  these  two 
N 
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different  methods.  The  values  denoted  with  ( ) are  those  caluclated 
with  given  bulk  Richardson  numbers,  and  those  without  { ) are  cal- 
culated based  on  given  h/L.  The  agreement  between  these  two  methods 
are  remarkable. 

The  geostrophic  drag  and  heat  transfer  coefficients  are 
defined  as,  (see  Yamada,  1976) 

Cjj  » k{(in(h/Z^)  - A)^  + 

= 1.35  k{£n(h/Z  ) - C}"' 

n o 

The  calculated  values  of  C,.  and  G based  on  the  two  different  methods 

D n 

# 

are  shown  in  Fig.  2.  Very  close  agreement  in  the  values  of  C and  C 

U n 

between  the  two  methods  are  indicated. 

B.  Surface  F 1 uxes  Ca 1 cu la  ted  by  the  Genera  1 i zed  Similarity 
and  Rossby  Number  Simi lari ty  Theor i es 

The  difference  between  the  generalized  similarity  and  the 
Rossby  number  similarity  theory  is  that  the  former  assumes  that  the 
boundary  layer  height  is  uniquely  determined  by  u,/f  and  L (Obukhov 
length).  In  the  generalized  version  h is  considered  as  an  independent 
variable  (Z i 1 i t inkev ich  and  Deardorff,  197^),  therefore,  in  the  latter 
the  effects  of  complicated  factors  such  as  nonstationari ty,  di'irnal 
heating,  large-scale  advectlon  of  heat  and  moisture,  large-scale 
subsidence,  etc.  on  which  h depends  in  the  real  atmosphere,  can  be 
considered  indirectly  by  specifying  h through  a rate  equation  (Deardorff, 

197^). 

In  this  section,  we  shall  examine  the  surface  fluxes  calculation 


based  on  these  two  theories.  For  this  purpose,  we  shall  use  the 
observed  data  set  taken  from  the  Wangara  experiment.  in  particular, 
those  days  of  the  Wangara  experiment  chosen  by  Melgarejo  and  Deardorff 
(1975)  will  be  selected  for  analyses.  The  observed  boundary  layer 
height  hg , defined  as  the  height  to  which  significant  cooling  had 
extended  as  judged  both  from  individual  profiles  and  their  evolution 
in  time,  will  be  used  to  compute  the  surface  fluxes  of  heat  and  momentum 
by  generalized  similarity  theory.  For  the  Rossby  number  similarity 
theory,  two  fixed  values  of  h will  be  examined,  i.e.  h * 1000  m and 
h = 500  m.  Determination  of  the  various  similarity  functions  are  done 
by  the  empirical  formulations  of  Yamada  (1976).  Specifically,  the 
similarity  functions  are  determined  by  the  bulk  Richardson  number 
approach  as  described  in  Section  3A. 

The  calculated  surface  friction  velocity  and  kinematic  heat  flux 

are  listed  in  Appendix  II.  The  Root  Mean  Square  Errors  (RMSE)  of  the 

calculated  surface  friction  velocity  u^.  and  kinematic  heat  flux  H /oC 

o p 

with  respect  to  the  values  calculated  in  Melgarejo  and  Deardorff  (1975) 
are  shown  in  Table  3.  It  should  be  noted  that  the  surface  friction 
velocity  and  kinematic  heat  flux  calculated  in  Melgarejo  and  Deardorff 
(1975)  are  based  on  the  flux-profile  relationship  of  Businger  et  al. 
(1971).  As  noted  in  Melgarejo  and  Deardorff  (1975),  the  fluxes  thus 
calculated  were  obtained  from  the  observed  wind  and  temperature  gradients 
between  1 and  k m.  The  fact  that  the  surface  similarity  theory  is  well 
established  and  the  profile-flux  relationships  of  Businger  et  al . (1971) 
are  well  deduced  from  a carefully  gathered  field  experiment  justifies 
the  use  of  the  results  of  Melgarejo  and  Deardorff  (1975)  as  the 
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comparison  standard.  From  Tabie  3 and  4,  we  see  that  during  the  unstable 
period,  the  RMSE  values  for  the  generalized  similarity  theory  are  not 
much  different  from  those  for  the  Rossby  similarity  theories.  This 
suggests  that  during  the  unstable  period,  the  Rossby  number  similarity 
approach  is  equally  valid  as  the  general  ized  similarity  theory.  During 
the  unstable  period,  wind  and  temperature  profiles  are  typically  well 
mixed  throughout  the  top  of  the  boundary  layer  (near  1500  m) . As  a 
consequence,  little  difference  in  the  prediction  by  the  three  similarity 
theories  may  be  expected.  During  the  stable  period,  it  is  evident  from 
Table  3 that  the  Rossby  number  similarity  theory  (especially  with  h = 

1000  m)  performs  better  than  the  generalized  similarity  theory.  This 
effect  may  be  due  to  the  fact  that  the  boundary  layer  heights  are  some- 
what indeterminate  at  the  stable  hours  so  that  the  computat ion  of  boundary 
layer  fluxes  based  on  the  generalized  similarity  theory  may  be  subject 
to  large  errors.  Based  on  this  analysis,  we  must  conclude  that  the 
Rossby  similarity  theory  is  useful  and  practical.  Although  the 
generalized  similarity  theory  has  more  of  a physical  basis  and  is  highly 
advocated  by  many  researchers  (see  e.g.  Deardor ff , 197^,  Arya,  1977, 
and  many  others),  the  advantages  of  the  theory  are  not  well  reflected 
in  these  surface  fluxes  calculation.  On  the  other  hand,  the  Rossby 
similarity  has  certain  practical  advantages  in  the  atmospheric  cir- 
culation models.  The  main  advantage  is  to  use  a fixed  scale  height  as 
the  top  of  the  boundary  layer  from  which  the  boundary  layer  fluxes  may 
be  computed. 

4.  SUMMARY  AND  CONCLUSIONS 

This  paper  concludes  that  the  height  of  the  constant  flux  layer  should 


not  be  assumed  to  be  greater  than  50  m in  a numerical  model.  This 
is  because  during  the  convective  periods,  the  very  super-adiabatic 
lapse  rate  normally  existing  near  the  ground  may  not  be  reflected  in 
the  calculation  of  a bulk  Richardson  number.  As  a result,  the  fluxes 
tend  to  be  underestimated.  On  the  other  hand,  if  the  height  of  a 
constant  flux  layer  is  assumed  to  be  very  shallow,  i.e.  less  than 
four  meters,  then  any  small  measurement  errors  in  the  wind  and 
temperature  profiles  will  lead  to  great  inaccuracies  in  the  calculation 
of  surface  shearing  stress  and  heat  fluxes.  For  numerical  models  of 
the  atmospheric  circulation,  we  recommend  a constant  flux  layer  of 
about  10-30  meters  depths  to  be  practical  and  reliable. 

We  describe  how  the  similarity  functions  A,  B,  C,  and  D can  be 
determined  based  on  the  bulk  Richardson  number  and  have  demonstrated 
that  the  similarity  functions  based  on  a recent  analyses  of  Yamada 
(1976)  can  be  computed  using  a Newton-Ral phson  iterative  scheme  during 
the  stable  period.  During  the  unstable  period,  a linear  relation  is 
found  between  the  bulk  Richardson  number  and  Obukhov  length.  The 
relationship  is  derived  and  the  results  of  drag  and  heat  transfer 
coefficients  are  found  to  be  satisfactory.  Analysis  based  on  the 
Wangara  data  indicates  that  the  Rossby  number  similarity  theory  is 
useful  and  practicable.  Although  the  generalized  similarity  theory 
has  a more  physical  basis,  this  analysis  shows  that  parameterization 
of  the  boundary  layer  fluxes  through  the  Rossby  number  theory  is 
equally  valid  as  is  the  generalized  theory  during  the  unstable 
period.  During  the  stable  period,  the  Rossby  number  similarity  theory 
performs  better  than  the  generalized  theory  in  the  computation  of  the 
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boundary  layer  fluxes.  However,  it  should  be  pointed  out  that  the 
Rossby  similarity  theory  may  not  be  valid  near  the  tropics  where 
the  determination  of  the  boundary  layer  height  by  u^^/f  is  no  longer 
valid.  For  this  reason,  the  generalized  similarity  theory  should  be 
recommended  for  modeling  the  global  circulations.  For  mid-latitude 
general  circulation  studies,  the  Rossby  number  similarity  theory  may 
be  adequate. 
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Table  1 : 


iFriction 
jVeloc  i ty 
(m  s"') 


1 


RMSE  (Root  Mean  Square  Errors)  of  the  Calculated  Surface 
Friction  Velocity  (m  s"M  and  Kinematic  Heat  Flux  (m  s"l  °K) 
for  the  Kansas  Experiment  Data 


Case  A (constant 

Case  B 

Case  C 

Case  D 

flux  layer  height. 

(CFL  -8m, 

(CFL  - 16  m. 

(CFL  - 32  m 

i . e. , CFL  “ 4 m; 

Z|  - 1 m. 

Z,  - 1 m. 

2,  - 1 m. 

and  Zj  « 1 m,  ■ 

4 m) 

Zj  * 8 m) 

- 16  m) 

- 32  m) 

0.054 

0.02k 

0.020 

0.018 

I 


Table  2: 


S imi lari 


ty  functions  A, 
h/Zo  = 10^ 


B and  C 


I 

1 


h/L 

A 

B 

C 

''iB 

(h/L) 

(A) 

(B) 

(c) 

-2 

1.90 

I.5A 

3.83 

-.12 

-2. A 

1.91 

1.A5 

3.87 

-10 

2.07 

0.93 

A. 38 

-.59 

-11. s 

2. 1 1 

0.88 

A.  50 

-18 

2.23 

0.77 

A. 82 

-1.03 

-21.0 

2.28 

0.73 

A. 95 

-26 

2.37 

0.68 

5.16 

-1.A6 

-29.7 

2.AA 

0.65 

5.30 

-38 

2.57 

0.60 

5.57 

-2.08 

-A2.A 

2.6A 

0.58 

5.70 

-60 

2.89 

0.52 

6.13 

-3.20 

-65.3 

2.96 

0.50 

6.2A 

-100 

3.35 

O.AA 

6.79 

-5.22 

-106.5 

3.AI 

0.A3 

6.88 

-li*0 

3.71 

0.39 

7.23 

-1.11 

-lAS.l 

3.78 

0.38 

7.31 

-180 

A. 00 

0.36 

7.55 

-9.33 

-190.2 

A. 07 

0.35 

7.62 

-250 

4.A1 

0.32 

7.96 

-13.01 

-265.2 

A.A8 

0.32 

8.03 

o 

A. 80 

0.29 

8.31 

-17. 8A 

-363.6 

A. 89 

0.28 

8.39 

-510 

5.32 

0.75 

8.75 

-27. 16 

-553.6 

5.A2 

0.25 

8.83 

Values  with  ( ) are  calculated  from  bulk  Richardson  number. 


le 


Table  3:  RMSE  (Root  Mean  Square  Errors)  of  the  Calculated  Surface  | 

Friction  Velocity  u*  (cm  s”')  and  Kinematic  Heat  Flux  :! 

Ho/pCp  (cm  s"'  °K)  with  Respect  to  the  Values  Calculated  i 

in  Melgarejo  and  Deardorff  (1975)  for  the  Wangara  Experiment  I 

Data. 


Generalized  Similarity 

Rossby  with  h » ,000  m 

Rossby  with  h « 500  ml 

U* 

Ho/dC 

P 

u* 

Ho/pCp 

Ho/pC  i 

P 

Unstable 

Period 

23.55 

6.3 

22.64 

4.5 

23.37 

1 

Stable 
iPer  iod 

, 

21., 



36.96 

9.8 

,4.6 

1 

19.52  1 

i 

r 


Table  Calculated  ^eans  and  Standard  Errors  of  Surface  Friction 

Velocity  u.  (cm  S“l)  and  Kinematic  Heat  Flux  Ho/cC 

I o ^ p 

(cm  s K)  for  the  Wangara  Experiment  Data 


Generalized  Similarity 

Rossby  with  h»1000  m 

1 

! 

1 Rossby  with  h“500 

Melgarejo  and 
Oeardorff  (1975) 

' 

Ho/oC 

P 

Ho/pC 

P 

i 

i u 

1 * 

B 

Ho/pC  1 

P ‘ 

Unstable 
!Per  iod 

214.  1 

27.95 

23.3 

23.89 

23.5 

28. 1 1 

25.2 

9.9  1 

, N-17 

(1.65) 

(6.70) 

(1.92) 

(5.28) 

(1.96) 

(6.55) 

(2.28) 

(1.30)  ! 

|stabl  e 
iPer  Iod 

2<4.0 

-2i*.98 

5.8 

-1  .6^4 

13.6 

-8.9 

10.2 

-.89  ' 

! N*6I 

1 

(2.62) 

(3.66) 



(1.63) 

(.73) 

(2.26) 

(2.32) 

(.73) 



(.08)  . 

1 

20 


j 


Bulk  Richardson  Number  (R.„) 

aB 

Fig.  2 — Geostrophic  drag  coefficient  Cq  (top)  and  geostrophic  heat  transfer  co- 
efficient Ch  (bottom)  calculated  by  given  h/L  values  (solid  lines)  and  by  the  Bulk 
Richardson  number  approach  (dashed  lines) 
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APPENDIX  II:  Surface  Friction  Veiocity  (cm  s *)  and  Kinematic 

' Heat  Flux  Ho/pC.  (cm  i ' °K)  Caicuiated  by  the  Generaiized 

P 

Similarity  and  Rossby  Number  Similarity  Theories  for  the 
Wan9ara  Experiment  Data 


Day 

Hour 

General ized 

Similarity 

Rossby, 

, h«l000  m 

Rossby, 

h»500  m 

Hetgarejo  & Oeardorff 
(1975) 

“* 

Ho/oC 

P 

u* 

Ho/pC 

P 

Ho/pC 

P 

Ho/cC 

0 

STABLE 

PERIOD 

1 

6 

14.5 

-7.49 

1.0 

-.06 

1.6 

-.16 

5.5 

-■35 

IS 

31.1 

-26.65 

6.5 

-.80 

14. 1 

-3.42 

11.0 

-1.24 

2i 

38.0 

-42.57 

11.8 

-2.62 

22.3 

-9.29 

12.4 

-1 .10 

24 

47.5 

-61.66 

8.6 

-1.75 

44.2 

-55-5 

16.5 

-1.42 

4 

3 

46.7 

-41.88 

11.4 

-2.42 

45.5 

-42.1 

14.2 

-1  .24 

6 

41.5 

-39.48 

10.0 

-1.97 

41.8 

-39.7 

15.8 

-1.77 

6 

18 

6.4 

-1.39 

0.2 

-.00 

0.8 

-.03 

5.8 

-.52 

21 

10.1 

-4.12 

0.4 

-.02 

1.4 

-.13 

3.3 

-.25 

7 

3 

7.5 

-2.95 

3.9 

-.71 

5.7 

-1.63 

3.9 

-.21 

6 

9.1 

-3.91 

2.1 

-.25 

3.6 

-■73 

5.0 

-.27 

iS 

5.6 

-1.12 

1.5 

-.09 

3.0 

-■31 

4.7 

-.55 

21 

3.2 

-.69 

2.7 

-.38 

2.3 

-■36 

1 1 .4 

-2.52 

24 

8.3 

-3.01 

2.3 

-.28 

3.9 

-■78 

10.4 

-1.36 

iO 

24 

6.4 

-.62 

0.3 

-.00 

0.7 

-.02 

0.8 

-.00 

ii 

24 

39.8 

-60.19 

1.7 

-.16 

7.7 

-1.94 

12.9 

-1.05 

12 

6 

10.9 

-4.28 

1.5 

-.14 

8.9 

-2.79 

15.4 

-1 . 18 

9 

33.7 

-26.39 

4.2 

-.53 

7.3 

-1.47 

21.6 

-.14 

21 

6.7 

-1.64 

3.4 

-.46 

6.1 

-1.38 

12.4 

-1.04 

24 

18.3 

-10.68 

4.6 

-.82 

7.9 

-2.29 

10.6 

-.76 

13 

3 

41.3 

-53.69 

2.7 

-.37 

50.0 

-04 . 80 

4.1 

-.23 

6 

50.3 

-105.80 

1.6 

-.17 

6.8 

-2.13 

9.8 

-.90 

IS 

31.7 

-30.96 

5.0 

-.61 

10.9 

-2.53 

4.4 

-.31 

21 

41.2 

-64.05 

5.7 

-1.03 

8.8 

-2.47 

1.4 

-.03 

24 

45.4 

-86.96 

2.8 

-.41 

12.8 

-5.54 

6.7 

-.61 

14 

3 

41.4 

-79.30 

1.4 

-.14 

8.1 

-2.70 

7.3 

-.66 

6 

41.7 

-67.10 

0.7 

-.05 

5.3 

-1.23 

13.1 

-1.49 

21 

13.4 

-6.05 

4.5 

-.73 

4.2 

-■77 

5.0 

-.43 

24 

16.3 

-9.19 

4.3 

-.80 

6.6 

-1.92 

8.0 

-.97 

i6 

24 

4.9 

-1.19 

3.1 

-.39 

1.6 

-■14 

10.6 

-1.18 

iS 

24 

37.9 

-38.07 

2.9 

-.29 

12.6 

-3^47 

13.5 

-1.25 

19 

3 

16.2 

-8.41 

2.3 

-.23 

4.4 

-■78 

11.4 

-■92 

6 

9.7 

-3.43 

1.6 

-.13 

2.5 

-■33 

10.0 

- . 68 

25 

21 

11.4 

-3.41 

0.9 

-.04 

3.0 

-■31 

3.4 

-.16 

26 

6 

4.2 

-.65 

1.4 

-.08 

1.6 

-.12 

3.8 

-.13 

21 

14.3 

-4.81 

11.6 

-2.17 

10.3 

-2  35 

6.4 

-.39 

30 

21 

45.6 

-39.00 

10.0 

-1.80 

48.0 

-42.00 

14. 1 

-1.26 

24 

50.4 

-53.11 

11.2 

-2.42 

47.0 

-53.27 

19.3 

-1.65 

4 
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APPENDIX  II  Cont'd 


Day 

Hour 

Generalized  Similarity 

Rossby, 

h«10Q0  m 

Rossby, 

h-500  m 

Melgarejo  £ Oeardorff 
(1975) 

Ho/oC 

P 

ho/pC 

P 

Mo/pC 

P 

Ho/pC 

0 

31 

3 

A5.9 

-55.30 

7.0 

-1.32 

17.9 

-6.79 

12.5 

-.92 

6 

AA.3 

-66.10 

A. A 

-.70 

21.8 

-10.86 

13. A 

-1.00 

21 

12.0 

-5.20 

3.5 

-.AS 

A. 7 

-.87 

9.0 

-.96 

2A 

32.8 

-58. 3A 

A. 3 

-.69 

7.5 

-2.01 

10.7 

-.83 

32 

3 

A.I 

-.88 

1.6 

-.15 

3.2 

-.55 

7.2 

-.AO 

6 

3.3 

-.81 

0.7 

-.OA 

2.7 

-.AS 

A. 6 

-.28 

2A 

10.0 

-A.  88 

0.8 

-.OA 

2. A 

-.33 

A. 3 
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